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Abstract 

We present tlie results of six Suzaku observations of the recurrent black hole transient 4U 1630 — 472 
during its decline from outburst from February 8 to March 23 in 2006. All observations show the typical 
high/soft state spectral shape in the 2-50 kcV band, roughly described by an optically thick disk spectrum 
in the soft energy band plus a weak power-law tail that becomes dominant only above ^^20 keV. The 
disk temperature decreases from 1.4 keV to 1.2 keV as the flux decreases by a factor 2, consistent with a 
constant radius as expected for disk-dominated spectra. All the observations reveal significant absorption 
lines from highly ionized (H-like and He-like) iron Ka at 7.0 keV and 6.7 keV. The brightest datasets also 
show significant but weaker absorption structures between 7.8 keV and 8.2 keV, which we identify as a 
blend of iron Kf3 and nickel Ka absorption lines. The energies of these absorption lines suggest a blue shift 
with an outflow velocity of ^ 1000 km s~^. The H-like iron Ka equivalent width remains approximately 
constant at ^ 30 eV over all the observations, while that of the He-like Ka line increases from 7 eV to 
20 eV. Thus the ionization state of the material decreases, as expected from the decline in flux. We fit the 
profiles with Voigt functions (curve of growth) to derive absorbing columns, then use these together with 
detailed photo-ionization calculations to derive physical parameters of the absorbing material. The data 
then constrain the velocity dispersion of the absorber to 200-2000 km s~^, and the size of the plasma as 
^ 10^° cm assuming a source distance of 10 kpc. 
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1. Introduction 

In recent years, a growing number of X-ray bina- 
ries have been found to exhibit absorption lines from 
highly ionized elements (e.g., Boirin et al. 2004 and ref- 
erences therein; Church et al. 2005). These systems 
range from microquasars such as GRO J1655— 40 (Ueda 
et al. 1998; Yamaoka et al. 2001; Miller et al. 2006) and 
GRS 19154-105 (Kotani et al. 2000; Lee et al. 2002) to 
low-mass X-ray binaries such as GX 13-1-1 (Ueda et al. 
2001; Sidoh et al. 2002), X 1658-298 (Sidoh et al. 2001) 
and X 1254-690 (Boirin & Parmar 2003). These systems 



are viewed at high inclination angles, and the absorp- 
tion structure is visible throughout the orbital period (e.g. 
Yamaoka et al. 2001; Sidoh et al. 2001; 2002). The absorp- 
tion features are therefore thought to originate in material 
that is associated with and extends above the outer accre- 
tion disk. This is illuminated by the X-rays produced from 
the innermost regions of the accretion flow (both the disk 
and hard X-ray coronal emission) . The reprocessed emis- 
sion and scattered flux from the extended wind can be seen 
directly in the Accretion Disk Corona sources, where the 
intrinsic X-rays are obscured (e.g. Kallman et al. 2003), 
but for the majority of highly inclined sources, the wind 
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material is seen in absorption against the much brighter 
intrinsic central X-ray source. Multiple absorption lines 
seen there then give an excellent probe of the physical 
conditions in the wind (e.g. Ueda et al. 2004; Miller 
et al. 2006), where the spectra indicate the presence of 
significant amounts of highly ionized outflowing material. 
Such accretion disk winds can be produced by thermal 
driving (e.g., Begelman et al. 1983), and/or by magnetic 
forces (Blandford & Payne 1982) as inferred by Miller et 
al. (2006) from Chandra HETGS data of GRO J1655-40. 

4U 1630—472 is a black hole candidate (e.g., Tanaka & 
Lewin 1995; McClintock & Rcmillard 2003), known for X- 
ray outbursts that repeat with an interval of roughly 600- 
690 days (Jones et al. 1976; Parmar, Angehni & White 
1995). The source has been observed with every major 
mission from EXOSAT (Kuulkers et al. 1997; Parmar et 
al. 1986), Ginga (Parmar et al. 1997), BATSE CGRO 
(Bloser et al. 1996), to Beppo-SAX (Oosterbroek et al. 
1998), RXTE (Tomsick et al. 1998; Kuulkers et al. 1998; 
HjcUming et al. 1999; Cui et al. 2000; Dieters et al. 2000; 
Trudolyubov et al. 2001; Kalemci et al. 2004; Abe et 
al. 2005), ASCA (Abe et al. 2005), and recently Chandra 
(J. Miller, private communication). These X-ray obser- 
vations showed that the source is absorbed heavily with 
TVh = (5-12) X 10^2 cm-2 (Tomsick et al. 1998; Kuulkers 
et al. 1997; Cui et al. 2000; Trudolyubov et al. 2001; Abe 
et al. 2005). While the IR counterpart was found during 
a 1998 outburst by Augusteijn et al. (2001), no optical 
counterpart is known for 4U 1630—472, most likely due 
to its high reddening and location in a crowded star field 
(Parmar, Stella & White 1986). Thus, the distance, incli- 
nation angle, and mass of the compact object are unknown 
for 4U 1630—472. This source is classified as a black hole 
not from dynamic mass measurements, but because of the 
similarity of its X-ray spectral and timing properties to 
those of systems with measured black hole masses. 

This source was observed six times with Suzaku during 
its most recent outburst as part of a program to study 
discrete spectral structures as a function of the changing 
accretion conditions. The Suzaku monitoring observation 
led to successful detection of absorption lines on all occa- 
sions. Though one of the absorption lines may have been 
detected by Chandra (J. Miller, private communication), 
the Suzaku XIS data reveal the line with a much higher 
significance, thanks to its larger effective area at 7 keV 
than that of Chandra. In this paper, we focus mainly on 
the absorption line structures, and how these discrete fea- 
tures evolve over the two months span of our observations. 
Considering the large absorption column density, we as- 
sume the source distance, D, at 10 kpc in this paper. In 
addition, by comparison to the other black hole binaries 
reported to show absorption lines, an inclination angle, 
I, of 70° is assumed. Errors quoted in this paper repre- 
sent 90% confidence limit for a single parameter unless 
otherwise specified. 



2. Observations and data reduction 

2.1. Observations 

An outburst of 4U 1630 — 472 was reported by the 
RXTE All Sky Monitor (ASM; Levine et al. 1996) on 2005 
December. From 2006 February 8 through March 23, we 
observed the source six times with a series of 21-23 ks 
exposures as part of the initial performance verification 
of Suzaku. Suzaku is the 5th X-ray astronomy satellite 
of Japan launched on July 10, 2005, with an M-V rocket 
from the Uchinoura space center (Mitsuda ct al. 2006). 
It carries 4 sets of X-ray telescopes (Serlcmitsos et al. 
2006) each with a focal-plane X-ray CCD camera (XIS, 
X-ray Imaging Spectrometer; Koyama et al. 2006) oper- 
ating in the energy range of 0.2-12 keV, together with a 
non-imaging Hard X-ray Detector (HXD, Takahashi et al. 
2006), which covers the 10-600 keV energy band with Si 
PIN photo-diodes and GSO scintillation counters. Three 
of the XIS (XISO, 2, 3) have front-illuminated (FI) CCDs, 
while XISl utilizes a back-illuminated (BI) CCD, achiev- 
ing an improved soft X-ray response but poorer hard X- 
ray sensitivity. Considering the large column density to 
4U 1630 — 472, we concentrate on the harder X-ray band in 
this paper, and only analyze the data from the FI CCDs. 

Figure 1 shows a lightcurvc of 4U 1630 — 472 obtained 
by the RXTE ASM. The Suzaku observations are indi- 
cated with downward arrows. Because the optical axes 
of the XIS and the HXD are slightly offset, the source 
was put at the center of the HXD field of view during the 
observations, to maximize the HXD count rate, and to 
reduce the XIS count rate as pileup is an issue for such 
a bright source (see below). The HXD was operated in 
the nominal mode throughout the observations. Because 
the source was very bright, the XIS was operated with 
the 1/4 window option, in which a smaller field of view of 
17'. 8 X 4'. 5 is read out every 2 s. A burst option was added 
in the first four observations. This decreases the effective 
exposure by a factor of two. The editing modes of the FI 
CCDs (XIS023) were set to 3 x 3 and 2 x 2 for high and 
medium data rates, respectively. Table 1 summarizes the 
details of the observations. 

2.2. Data reduction 

For both the XIS and the HXD, we use version 0.7 
screened data provided by the Suzaku team. The screen- 
ing of the version 0.7 XIS data are based on the following 
standard criteria: a) Only GRADEO, 2, 3, 4, 6 events are 
accumulated, b) The time interval after passage through 
the South Atlantic Anomaly is larger than 436 seconds, 
c) The object is at least 5° and 20° above the rim of the 
Earth during night and day, respectively. We limit our 
analysis to the 3x3 mode data (high data rate) , because 
the 2x2 mode data (medium data rate) are not well cali- 
brated at this time. During the 2nd and 5th observations, 
more than half of the exposure was acquired in 2 x 2 mode, 
and thus these data have poorer statistics. The exposure 
times for the total and high data rate are summarized in 
Table 1. 

The XIS events were extracted from a circular region 
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Epoch Start End XIS options Exposure [ks] * Flux [erg s ^ cm •^J 

totalt Hight 2-10 keV 20-40 keV 



1 


Feb 8 8:14 


Feb 9 3:20 


1/4 window. 1-s 


burst 


22.2 


22.2 


4.6 X 10' 


-y 


6.9 X 10' 


'ii 


2 


Feb 15 23:08 


Feb 16 13:08 


1/4 window, 1-s 


burst 


21.5 


9.6 


4.0 X 10" 


-9 


6.7 X 10' 


-11 


3 


Feb 28 23:00 


Mar 1 16:43 


1/4 window, 1-s 


burst 


21.3 


21.3 


3.1 X 10" 


-9 


5.0 X 10' 


-11 


4 


Mar 8 01:39 


Mar 8 17:31 


1/4 window, 1-s 


burst 


21.2 


20.7 


2.8 X 10' 


-9 


3.5 X 10- 


-11 


5 


Mar 15 5:19 


Mar 15 17:51 


1/4 window 




23.2 


5.2 


2.5 X 10' 


-9 


2.6 X 10" 


-11 


6 


Mar 23 9:57 


Mar 23 22:21 


1 /4 window 




21.7 


15.4 


2.2 X 10' 


-9 


1.5 X 10' 


-11 



* Without correction tor the burst option. With the burst option, the net exposure is a halt ot the nominal exposure listed 
in this table. 

t Not exposures of total (medium and high data rates) and only high data rate are shown. For FI-CCDs, observations were 
performed with 3x3 and 2x2 mode in data-rate High and Med, respectively. 
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Fig. 1. The 1.5-12 keV light curve of 4U 1630 -472 (top 
panel), and hardness ratio of 5-12 keV to 3—5 keV (bottom 
panel), obtained with the RXTE ASM since 2005 August 13 
(MJD 53595). The Suzaku observations are indicated with 
down-arrows. 

with a radius of 4'. 3 centered on the image peak. This 
extraction circle is larger than the window size. The ef- 
fective extraction region is therefore the intersection of the 
window and this circle. The background was entirely neg- 
ligible throughout the six observations (less than 1% of 
the signal at 9 keV), so was not subtracted from the data. 
Using XIS2 as a reference, we find that the 2-10 keV count 
rate declined over the six observations, with 124, 112, 90, 
82, 74, and 64 cts respectively, roughly corresponding 
to 4.6, 4.0, 3.1, 2.8, 2.5, and 2.2 x 10"^ erg s"! cm'^ in 
the same energy band. Within each observation, we find 
no significant intensity variation intrinsic to the source. 

The source was so bright during these observations that 
the XIS suffered from photon pileup at the image center. 
A rough threshold above which the photon pileup becomes 
significant is ~100 cts exposure"^ for a point source. The 
count rates of 4U 1630 — 472 during the observations were 
around this level. Furthermore, because 4U 1630 — 472 has 



a very soft energy spectrum, the effect of photon pileup 
can be significant in the hard band spectra. Therefore, 
we excluded a circular region of 30" at the image center 
from the event extraction region to minimize these effects. 
This did not remove the photon pileup completely, and an 
artificial structure is seen as a hard excess above ^ 8 kcV 
at the <5-10% level. Because photon pileup tends to 
produce X-ray events which have twice (or a little more) 
as high an energy as the most abundant events (in this 
case, about 4 keV), we modeled the artificial structure as 
a gaussian centered at 10 keV with a fixed width of cr = 
1 keV. As seen in the next section, this phcnomcnological 
model works quite well to compensate for the effect of 
pileup. 

Because the event extraction region has a complex 
shape (an annulus truncated with the window size), 
we calculated the effective area for each XIS sensor by 
using XISSIMARFGEN (version 2006-04-24), which calcu- 
lates the arf file through Monte Carlo simulations. The 
three FI CCDs have almost the same performance, so 
we have added their data to increase the statistics in 
the spectra. Hereafter wc refer to the summed data as 
XIS023. We calculated the corresponding instrumental 
responses by summing the latest redistribution matrices 
of ae_xi[0,2,3]_20060213.rmf using the ftool ADDRMF and 
the ARF files using the ftool ADDARF. 

Though the main objective of this paper is the absorp- 
tion structure in the XIS energy band, it is also useful to 
understand the continuum emission since this provides a 
knowledge of spectral states and photoionising flux. We 
thus analyzed the HXD data using only the PIN data, 
since the source is dim in the GSO band. Both high and 
medium data rate are used in the PIN data analysis. The 
version 0.7 PIN data are screened by the following se- 
lection criteria: a) The object is at least 5° above the 
earth rim, b) The time interval after passage through the 
South Atlantic Anomaly is longer than 500 seconds, c) 
The cutoff rigidity is greater than 8 GeV c~^. The PIN 
background events were provided by the HXD team for 
each observation based on two different models, A and 
B (Kokubun et al. 2006). Here we constructed the PIN 
background spectra based on model B filtered by the same 
good time intervals. The latest version of the response 
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matrix, ae_lixd_pinhxnom_20060419.rmf, was used in our 
spectral analysis. Even if we allow 5% systematic uncer- 
tainty in the background level, the source is clearly de- 
tected up to 40 keV and 35 keV in PIN data from the 
first four and last two observations, respectively. Dead 
time fractions were calculated from the pseudo event rates 
(Takahashi et al. 2006; Kokubun et al. 2006) as 5.3, 4.3, 
4.6, 5.1, 4.6, and 4.6%, in the first to last observations, 
respectively. After the dead time correction, the back- 
ground subtracted PIN count rates are estimated to be 
0.16, 0.13, 0.093, 0.065, 0.055, and 0.043 cts s'^ in the 
range of 20-40 keV, corresponding to 6.9, 6.7, 5.0, 3.5, 
2.6, and 1.5 x 10^^^ erg s~^ cm~^. 

2.3. Calibration uncertainties of the XIS data 

Because the XIS calibration is still in progress, we in- 
clude a brief description of the calibration uncertainties. 
The XIS energy scale is relatively well calibrated with an 
uncertainty of less than 5 eV, or 0.2% when the standard 
mode is used (Koyama et al. 2006). However, when a 
window option is applied, the gain can change slightly, 
because the number of frame-store and parallel transfers 
are different than in the standard mode. This difference 
is not taken into account in the version 0.7 processing 
software. Since a window option was used for all the ob- 
servations of 4U 1630 — 472, the XIS instrument team es- 
timates that this produces a systematic gain difference 
of ~10 eV at 6 keV at the time of 4U 1630 - 472 obser- 
vations. When the window option is used, the onboard 
^^Fe calibration source is not directly visible but some 
events can be scattered into the active window in XIS2, 
which has the strongest calibration sources. The scat- 
tered ^^Fe events are accumulated from the regions free 
from the source photons (1/8 of the window area for each 
side). We obtain an average energy of 5891 ± 20 eV. This 
is consistent with the ~10 eV offset estimated above. We 
also examined the data for a possible gain shift due to the 
large surface brightness (local event density) of the image, 
by comparing the absorption line energies in the observed 
spectra for different radii of extraction annuli. We find 
no significant differences. Thus a conservative estimate is 
that the XIS gain at the iron band is correct to within 
20 eV during the 4U 1630 - 472 observations. 

The energy resolution of the XIS is gradually degrad- 
ing due to the radiation damage on the CCD. However, 
this effect is not yet included in the currently available 
response matrices. At the time of the 4U 1630 — 472 ob- 
servations, a narrow line at 6 keV would appear broad- 
ened with cr ~ 40 eV, where a is an equivalent gaussian 
width (standard deviation). Thus, when we evaluate a 
line width, we take this effect into account. It is impor- 
tant to note that this does not affect the measurement of 
the equivalent width. 

The effective areas of the XIS and PIN are calibrated us- 
ing observations of the Crab Nebula. The accuracy is ex- 
pected to be 5-10%. This is sufficient to estimate various 
emission parameters of 4U 1630— 472. However, because 
of the very good statistics of the data, the systematic er- 
rors in calibration could degrade the model fitting signif- 



icantly. A major contribution to the systematic error is 
the calibration uncertainty at gold M-II edge in the instru- 
mental response. To account for this residual, we added a 
gaussian at ~ 3.2 keV with a fixed width of cr = 0.1 keV. 

3. Analyses and results 

3.1. Broad band features 

Figure 2a shows the 2-50 keV spectra of 4U 1630 - 472 
obtained with the XIS023 and the PIN. We do not use 
data below 2 keV, because of large interstellar absorp- 
tion. An instrument-independent way to view the spectra 
is to construct the ratio of the source spectra to the cor- 
responding data from the Crab Nebula, which has an ap- 
proximately power-law spectrum with F = 2.1, and inter- 
stellar absorption of A^h = 3 x 10^^ cm~^ (Toor & Seward 
1974; Schattenburg & Canizares 1986). Figure 2b shows 
the Crab ratios calculated using the Crab data obtained 
with Suzaku on 2006 April 5 at the HXD nominal posi- 
tion and the XIS operated with the 1/4 window option. 
This plot reveals several noticeable features in the broad 
band source spectra. The source spectra are characterized 
by a dominant soft thermal component in the XIS band, 
which is heavily absorbed, while the harder X-ray PIN 
data clearly show a weak power-law tail. 

The dominant soft component is generally observed 
from bright black hole binaries in the high/soft state 
(Tanaka & Lewin 1995), and is believed to be emission 
from the optically thick standard accretion disk (Shakura 
& Sunyaev 1973). In the Crab ratio (figure 2b), the weak 
power-law tail is roughly characterized as F ~ 2 in the 
first two observations. When the 20-50 keV PIN spec- 
tra of these two observations are fitted with the single 
power-law model under the current response matrix, the 
values of F are obtained as 2.4 ±0.4 and 2.5 ±0.5, respec- 
tively. The extrapolated 20-100 keV fluxes are estimated 
as (1.4±0.3) X 10-^° and (1.1±0.3) x 10"^'' erg s'^ cm'^, 
respectively, which are only 3 % of the absorbed disk fiux 
in the range of 2-10 keV. Though the steeper power-law 
tails may be suggested on the latter four observations 
from the Crab ratio plot (figure 2b), they are difficult 
to quantify because of their faintness under the current 
background uncertainty. Further detailed analysis of the 
HXD data will be presented in a separate paper. 

The dominance of the disk and the weak power-law tail 
are both characteristic of the high/soft state. In addition, 
the Crab ratios also reveal the presence of complex ab- 
sorption structures in the iron K band. In this paper, we 
concentrate on these absorption line structures and their 
relation to the accretion disk parameters. Hereafter, we 
perform detailed spectral analyses on the XIS data. 

3.2. XIS spectra with canonical multi-color disk and neg- 
ative gaussian fits 

Figure 3 shows the 2-10 keV XIS raw spectra. To repro- 
duce the continuum spectra, we use the multi-color disk 
model (diskbb in xspec\ Mitsuda et al. 1984; Makishima 

^ http:/ /hcasarc. gsfc.nasa.gov/docs/xanadu/ 
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Fig. 2. (a) The X-ray spectra of 4U 1630-472 obtained with the XIS023 and the HXD PIN. The first to the last observations are 
shown with different colors as black, red, green, blue, right blue, and magenta, (b) The ratio of 4U 1630 — 472 spectra to the Crab 
spectrum. 



et al. 1986) modified by interstellar absorption (wABS 
model; Morrison & McCammon 1983), following the 
canonical modeling for the high/soft state black hole bina- 
ries. Here, the diskbb model approximates emission from 
the standard accretion disk, and is known to successfully 
reproduce the dominant soft spectra of the high/soft state 
black holes. As shown in figure 2b, the contribution of the 
power-law tail is negligible in the XIS energy band, so is 
not included in the spectral fitting. The residual effects 
of pileup are modeled by a broad gaussian at 10 keV (see 
§ 2.2), and the current uncertainties in the instrument 
response at the gold M-II edge by another gaussian at 
3.2 keV (see § 2.3). 

We start with the first observation, which showed the 
brightest spectrum. The continuum is well reproduced 
by the diskbb model with a disk temperature of kTia ~ 
1.39 keV. This is amongst the highest disk tempera- 
tures observed with RXTE from this source in previous 
high/soft states (e.g., Abe et al. 2005). A disk bolometric 
luminosity of Ldisk = 2.8 x 10'^* • a erg s~^ and an apparent 
disk inner radius of rjn = 24 • ^/a km are obtained, where a 
is defined as {D/IQ kpc)^/(cosj/cos70°). The absorbing 
column density, A^h, is estimated to be 8.15 x 10^^ cm~^. 
This is close to that seen with previous observations by 
ASCA and RXTE (e.g., Abe et al. 2005). 

The fit with this simple model is however far from ac- 
ceptable {x^ /dof ~ 1465.0/574) because of the complex 
structure in the iron K band. Residuals between the data 
and the best-fit continuum model are shown in the mid- 
dle panel of figure 3a. The most obvious features are 
two narrow dips at 7.0 keV and 6.8 keV. We therefore 
add two negative gaussians to the model. The fit is dra- 
matically improved with x^/dof = 714.0/570. The line 
center energies are estimated to be 7.00llg QQ5 keV and 
6.73±0.02 keV, consistent with the H-like and He-like iron 
Ka lines at 6.966 keV and 6.697 keV ^; the lines are possi- 



bly blue-shifted with 2 ~ 5 x 10"'^. The equivalent widths 
are 30 eV and 7 eV, for the H-like and He-like lines, re- 
spectively. The lines are narrow with an upper limit to <j 
of 40 eV and 50 eV. These values are consistent with what 
we expect for narrow lines given the current response ma- 
trix, which does not include the degradation of the energy 
resolution (§ 2.3). We thus fix the value of a at 10 eV to 
represent intrinsically narrow lines in this subsection. 

In addition to these obvious line structures, the resid- 
uals reveal weaker absorption structures at 7.8 keV and 
8.2 keV. These energies are very similar to the rest frame 
energies of He-like and H-like iron K/3 lines, namely 
7.88 keV and 8.25 keV. However, they are also close 
in energy to the He-like and H-like nickel Ka lines at 
7.80 keV and 8.09 keV, respectively, which can have com- 
parable equivalent width to that of the iron K/3 lines 
(Kotani et al 2000). The data can marginally constrain 
two weak lines at these energies, but not four. Hence we 
include two additional negative gaussians, again with the 
width a fixed at 10 eV, and allow the energies to vary 
freely. This decreases y^/dof to 686.4/566, indicating 
that the two negative gaussians are statistically signifi- 
cant at the 99.98% confidence level using an F-test (where 
Fiiy^Av)- ( Ax^ /Aiy) /xl = 5M). The best-fit parameters 
are summarized in Table 2. Residuals between the data 
and the best-fit model are shown in the bottom panel in 
figure 3a. The inferred energies are almost consistent with 
H-like and He-like iron K/3, though they can be contami- 
nated by nickel Ka. 

We performed the same spectral analysis on the other 
five observations. The XIS spectra, the best-fit models, 
and the residuals are shown in figure 3b-f. The best- 
fit continuum and absorption parameters are summarized 
in Table 2, and the time histories of the continuum pa- 
rameters are given in figure 4. Through the six observa- 
tions over two months, the value of rin remained approx- 



Weighted averages of K«i and Ka2 by referring to the 
National Institute of Standards and Technology (NIST; 



see http://aeldata.phy.nist.gov/archive/data.html ) and Drake 
(1988) for H-like and He-like Ka, respectively. 
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imately constant at 25 • ^fa km, while the disk tempera- 
ture gradually decreased to 1.18 keV, corresponding to a 
factor 1.75 decrease in the disk bolomctric luminosity to 
idisk = 1-6 X 10"^* • a erg s~^. The constancy of is an- 
other characteristic feature of the high/soft state of black 
hole binaries (e.g., Tanaka & Lewin 1995 and references 
therein), and confirms that the object was indeed in the 
high/soft state. The correlated decrease in the derived 
value of iVn by ~3% could be an artifact of the over- 
simplified DISKBB model; it assumes that the disk radiates 
locally as a blackbody, but distortions from this are sug- 
gested both observationally (e.g., Kubota & Makishima 
2004; Abe et al. 2005) and theoretically (Watarai et al. 
2000; Davis et al. 2005) from such luminous binaries. 

Irrespective of the continuum modeling, the two strong 
iron Ka lines are clearly detected throughout the obser- 
vations, whereas the weaker lines at higher energies are 
significantly detected in the first four datasets (F-test sig- 
nificance values are given in Table 2). Figure 5 shows an 
enlargement of the 6-9 keV spectra, and residuals between 
the data and the continuum model. The time histories 
of the individual absorption-line equivalent widths and 
their energies are plotted in figure 6. While the equiva- 
lent width of the H-like line was almost constant through- 
out the observations, that of the He-like line increased 
by a factor of two between the first and second observa- 
tions, and then perhaps increased slightly more through 
the subsequent observations. These imply that the ion- 
ization state of the absorbing gas decreased over time, as 
expected from the declining source luminosity. 

The time averaged center energies of the two strong 
iron Ka lines are estimated to be 6.987 ± 0.005 keV 
and 6.714 ± 0.009 keV, suggesting a blue shift veloc- 
ity of 900 ± 200 km s"! and 700 ±400 km s"!, respec- 
tively. If we take only the first observation, the blue 
shifts become more significant with 1500^2oo ^^'^ 
1500 ±900 km s~^, for the H-like and He-like irons, re- 
spectively. Here, the errors are statistical only. Given the 
current uncertainty in the gain calibration of ^ 20 eV or 
~ 900 km s~^ (see § 2.3), the blue shift cannot be mea- 
sured with higher significance. However, a blue shift is 
physically expected, since the material responsible for the 
absorption lines is certainly a wind from the disk (§ 1), 
and these values are consistent with the observed blue 
shift from other accreting systems like GRO J1655 — 40 
(Miller et al. 2006) and GX 13 + 1 (Ueda et al. 2004). In 
the next subsection, the absorption line parameters and 
their time evolution are studied in detail. 

3.3. Reanalysis of the absolution lines with Voigt profile 

The simple gaussian analysis above shows significant 
He-like and H-like iron Ka absorption lines in the data, 
and suggests the presence of the weaker absorption struc- 
tures at 7.8-8.2 keV. Since the parameters of the weaker 
high-energy lines are difficult to constrain, below we fo- 
cus on the stronger two lines, and study them using 6.0- 
7.5 keV spectra. The limited energy range means that we 
do not need to include the (instrumental) broad gaussian 
at 10 keV or the narrow gaussian at 3.2 keV. We model 
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Fig. 4. Time histories of continuum parameters. Values of 
^diski kTin, Hn j and x^/d.o.f are plotted. As described 
in the text, the decrease of Nj^ could be caused artificially by 
the over-simplifiod DISKBB model. 

the continuum by DISKBB as before, with free parameters 
Tin and Hn, but fix TVh to 8.0 x 10^^ cm~^, which is the 
average value obtained from the 2-10 keV fits. This does 
not affect the derived line center energies or equivalent 
widths, since the differences of A'h from the fixed value, 
at most ^ 1 X 10^^ cm^^, produce a rather small change 
to the iron edge at 7.1 keV. 

To conduct a detailed analysis of these absorption lines, 
we use a full Voigt profile as implemented in the kabs 
model. The kabs model is constructed by Ueda et al. 
(2004) (see their §3.2.1) as a local model for xspec, and 
has three fitting parameters for each specific ion; the ion 
column density, A^ion, the velocity dispersion, 6, and the 
bulk outfiow velocity characterized by a blue shift, z. 
The line shape is produced through a convolution of a 
Lorentzian profile, describing natural and pressure broad- 
ening with a gaussian describing Doppler motions (spec- 
ified by the velocity dispersion including thermal and ki- 
netic motion). Except for forbidden lines, each individual 
line of the Ka multiplet (Kai and Ka2) is included sep- 
arately, with the appropriate energy, oscillator strength, 
and Einstein coefficient. This properly models the fine 
structure in the line even though these details of the line 
profile cannot be resolved by the CCD. The values of chi- 
squared are the same as for the previous negative gaus- 
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Fig. 3. XIS023 spectra for the first to the last observations in panel (a) to (f) in this order. The observed data and the best-fit 
DISKBB model with negative gaussians are shown in the top panels. The middle panels show residuals between the data and the 
model without negative gaussians. The bottom panels show the residuals between the data and the best-fit model. Note that the 
burst option was not used in the 5th and 6th observations (Table 1), hence the pileup effect is more obvious in panels (e) and (f) 
in spite of the gradual decrease of the source fiux. An edge-like structure may be indicated at ~ 4 keV of which residual level is 
within 1%. This structure is understood as a coupling of calcium K edge at 4.05 keV of interstellar absorption and the uncertainty 
of current response matrix. 
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Table 2. Best-fit parameters for the 2-10 kcV XIS spectra 



epoch 


continuum 


gaussian 1 


gaussian 2 


gaussian 3 


gaussian 4 


x'/dof 




A^H [10^2 cm-2] 
kTi^ [kcV] 


£^c [keV] 
EW [cV] 


Ec [keV] 
EW [eV] 


Ec [kcV] 
EW [cV] 


Ec [keV] 
EW [cV] 






Tin [\/a km] 


[a [cV]) * 


{a [eV]) * 


significance under i^-test ^ 




1 


8 

°-^3-0.03 
1 '}er;+0.004 

1.00d_Q QQ3 


6.73 ±0.02 

'-2 


/.UUl_o.005 

30 ±1 


7.87 ±0.04 
6.6^^Coi)* 


8.24±0.05 


686.3/566 




24 


(<37) 


(<49) 


F(566,4) = 5.69 


99.98% 




2 


8.09 ±0.04 


' ^-0.02 

17±3 


6 005+0.009 
34ti 


7.80 ±0.05 
11.4±5(l°i) 


23.1^r(tg5) 


628.6/566 




24.2l0| 


(<46) 


(13111) 


F(566,4) = 9.68 


> 99.999% 




3 


8.06 ±0.03 


6.72 ±0.01 


u.»0'±_o 006 


7.77 ±0.08 


8.23 ±0.04 


538.8/566 




1.262 ±0.003 


18 ±2 


32 ±2 


7.74±4(t%4) 








rjr 9 + 0.2 
^"-'■^-0.2 


(<37) 


(18±18) 


F(566,4) = 6.89 


99.998% 




4 


7.97 ±0.03 

i.Z4U_Q 004 


6.72 ±0.01 
21 ±3 




7.84±1!:- 
18.2t^(tS.i) 


-, -, +U.Ob 
15.8_7(_o.2) 


636.7/566 




24.9l°J 




(<66) 


F(566,4) = 12.0 


> 99.999% 




5 


7.94 ±0.04 


6.72 ±0.02 


6.97±0.01 


(7.8-7.9)?! 


(8.2-8.3)?* 


599.9/566 




1 99n+"0"4 
i.zzu_o.o05 


20 ±4 


28+5 


< 13 


< 10 






24-9lo.3 


(<55) 




F(566,4) = 0.34 


15% 




6 


7.92 ±0.03 


6.705 ±0.009 


fi qfiq+0.008 
o.yDy_o.oo6 


7.85 ±0.09 


(8.2-8. 3)§ 


616.0/566 




1.183± 0.003 


22+^ 


31 ±2 




<9 






1+0.2 
^"J-^-0.2 




(<23) 


i^(566,4) = 1.05 


62% 





* The best-fit values of cr, or the upper Umit if the best-fit value is smaller than 10 cV. 
t F-values are presented for the addition of gaussians 3 and 4. 

•I-The best fit normalization factors of the broad gaussian at 10 keV for the pileup are 7.0, 4.1, 2.4, 1.9, 5.2, and 3.1 X 
10~^ photons s~^ cm~^, in order of the first to the last observations. Uncertainty of the pileup effect to the equivalent 
width is considered by including/excluding the broad gaussian to the corresponding continuum. 

§ Line center energies are fixed at those of He-like and H-like iron K/3 in the rest frame to obtain upper limits of the 
equivalent width, since they are not constrained by the fitting. 



sian fits, but this model gives the ion column A^ion tor 
an assumed velocity dispersion, b, as opposed to simply 
an equivalent width. We require the two lines to have 
the same b and the same z, and neglect any line emission 
which could fill in the absorption. 

Although we cannot measure b from line profiles, its 
upper and lower bounds may be obtained. The upper limit 
is set at 6~ 2000 km s""'^ by the upper limits on the line 
width, fT ~ 30-40 eV (§ 3.2). The lower limit to b is given 
by the thermal velocity. At such high ionization, the self- 
consistent temperature of the photoionized material must 
be of order the Compton temperature, Tq. Considering 
the balance between Compton heating by the dominant 
disk illumination, and Compton, bremsstrahlung and ion 
coolings, Tq is roughly estimated as /cTin/3 (i.e., 0.5 keV 
in this case). The flux in the power-law tail has only a 
very small effect on this as it is so much weaker than the 
disk (20-100 keV power-law flux is only ^ 1% of the disk 
bolometric flux), and docs not contribute significantly on 
the compton heating (e.g.. Reynolds & Wilms 2000). This 
temperature implies thermal (random kinetic) velocties 
as 6 = y^2kTc/ma,tom, and thus = 42 km s^^ can be 
obtained as a secure lower limit. 

Table 3 gives the derived Nion for these two limiting 



velocities, together with velocities of 5 = 100, 200, and 
500 km s^^. As expected from the curve-of-growth anal- 
ysis, the observed equivalent width can be produced by a 
lower value of 7Vio„ if a higher value of b is assumed. In 
Table 3, we can put further constraints on b by using the 
observed equivalent widths of weaker lines at 7.8 kcV and 
8.2 keV (gaussian 3 and 4 in Table2). The table includes 
the model predictions for the equivalent width of the iron 
K/? absorption lines. There, models with "1" in the last 
column ("Note") are inconsistent, because the predicted 
iron K/3 line equivalent widths exceeds the actually ob- 
served values, leaving no room for the nickel Ka lines to 
contribute. By scanning b with finer steps, wc found that 
6 of ^ 60-200 km s~^ significantly over predicts the K/3 
equivalent width, and is hence ruled out. 

The top three panels of figure 7 show the time histories 
of the best-fit kabs parameters, assuming b = 500 km s~^ 
which is based on the marginally resolved line width in the 
Chandra HETGS data of GX 13 + 1 (Ueda et al. 2004). 
The column density of the H-likc iron, N-pa xxvi: is almost 
constant at ~ 1 x 10^* cm~^ while that of the He-like iron, 
A'fc XXV I increases significantly from the first to the sec- 
ond observation. This is consistent with the behavior of 
the equivalent width estimated using the negative gaus- 
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Fig. 5. The XIS023 6-9 keV spectra concentrating on the ab- 
sorption hne structure. Data and the best-fit models (Table 2) 
are shown in the top panel. Residuals between the data and 
the best-fit model excluding the four negative gaussians are 
shown in the successive six panels. Colors are the same as in 
figure 2. 

sian fit (Table 2). The weighted average blue shift is about 
z = (3.3±0.7) X 10~^, corresponding to an outflow velocity 
of 1000 ±200 km s~^. Again, the first observation shows 



(5.0 



+0.8 



0.6)X10" 



or V ■- 



1500 



-1-240 
-180 



km s ^ (see Table 3). 



They are consistent with the results obtained with the 
negative gaussian fits in § 3.2. 

4. Physical parameters of the absorber 

We observed very significant absorption lines from He- 
like and H-like iron Ka in the disk dominated X-ray spec- 
trum of 4U 1630 — 472 during its decline from outburst. 
The inferred blue shift of the absorption lines supports 
the presence of an outflow. This source thus joins the 
growing number of Galactic binaries with such absorp- 
tion features, suggesting that winds are a generic feature 
of bright accretion disks observed at high inclination an- 
gles. Through the analyses in § 3, we found that iron is 
always predominantly H-like, requiring that the ioniza- 
tion parameter, ^ = L/nr^ (where n and r are the number 
density and the distance of the absorbing material from 
the illuminating source), be very high. In this section. 
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Fig. 6. Time histories of absorption line parameters based on 
the negative gaussian fit in the range of 6—10 keV. The first 
two panels show line center energies and equivalent width of 
gaussian 1 (filled circle) and gaussian 2 (open circle). The 
bottom two panels show those of gaussian 3 (filled triangle) 
and gaussian 4 (open triangle). On the top panel, average 
values arc shown with dashed lines. 

we use the xstar photoionization code'^ (version 2.1kn5, 
Kallman & Bautista 2001) to calculate the ionization bal- 
ance of the line-producing material, to convert the ion 
column densities into physical parameters including total 
column density, number density, and location, for each 
observation. 

4-.1. First observation 

We derive the physical parameters of the absorber using 
the first observation as a specific example to illustrate 
the method. We calculate the relative ion populations 
under illumination by a DISKBB spectrum, and use this to 
predict a theoretical ratio of TVpc xxv to iVpo xxvi- The 
observed ratio from the data then enables us to estimate 
the ^-parameter, and so convert the observed ion column 
densities of H-like and He-like irons into a total column 
density, iVtot- 

We use the illuminating spectrum derived from the 
Suzaku data, i.e., a diskbb spectrum of kTi^ = 1.4 keV, 
to calculate the self-consistent ion populations of iron us- 
ing XSTAR. The solid lines in figure 8a shows the cal- 
culated relative fractions of each ionization stage of iron, 
/ion, against ^, while figure 8b shows Npc xxv/^Fc xxvi- 

http:/ /heasarc. gsfc.nasa.gov/docs/softwarc/xstar/xstar. html 
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Table 3. Best-fit parameters of tlie 6—7.5 keV KABS model 



epoch 



The best-fit KABS parameters 



.^Fe XXV 

IQis cm"2 



iVpc 



XXVI 



10-3 



predicted K/3 EW* 
Fe XXV Fe XXVI 

(96)# cV eV 



A^ion ratioT 
He/H 
xlO-2 



plasma parametcrs-f Note 
10"^ 10^3 cm"2 



& = 42 km s-^ (fcTkin = 0.5 keV) 



1 


^•"-0.3 


34+5 


O-^-0.7 


116.1 


4.1 


12.5 


i.7ti:S 


19+1'^ 


-,on+i4() 


2 


2 




40+1° 


4.2ti-; 


97.9 


9.3 


11.3 




5tt 


on+39 


2 


3 






o + l.O 


98.8 


7.9 


12.6 


14+9 


64 




2 


4 


6.7;^;3 


3ii^ 


2 4+1-3 


102.4 


10.1 


11.8 


224' 


5t? 


22+12 


2 


5 


5.7t3j 




1.5±1.5 


85.7 


8.1 


10.3 


01+28 
^^-12 


5+5 

3-2 




2 


6 


7 7+2.5 
'■'-1.8 


33+4 


1 q+i-o 


108.2 


10.3 


12.0 


24ti'5 


6±2 


22t^ 


1,2 



b = 100 km 



^ (fcTkin = 2.9 keV) 



w 

4 

+ 10 

-9 

23t^ 

+7 
-6 
+ 10 

8 
+7 
-6 



5 3+'^- '^ 


115.9 


3.1 


17.1 


1 0+1'^ 


294i 


140+™ 


1, 2 


4.4+1-; 


98.0 


11.3 


17.6 


7+9 
' -4 


84 


45+^^ 


2 


o n+0.9 
-^■^-0.8 


99.2 


11.1 


17.1 


91^ 


q+6 


26t?° 


2 


9 7+1.2 
^- ' -0.8 


103.5 


12.9 


16.8 


154' 


6+^ 

"-3 


17+17 


2 




87.5 


11.9 


16.4 




7+8 
' -4 


14+f 


1, 2 


9 q+0.8 
^•■-•-0.9 


107.2 


13.9 


17.0 


184' 


7+4 
' -2 


16t- 


1, 2 



1 


0.23 


2 


2.03 


3 


1.98 


4 


3.12 


5 


2.51 


6 


3.99 



-0.09 
+ 1.30 
-0.80 
+0.97 
-0.78 
+ 1.46 
-1.07 
+2.17 
-1.17 
+ 1.62 
-1.21 



24 
31 



21 
17 
22 



b = 200 km s-^ (fcrkin = 12 keV) 



1 


^•J^^-0.04 


3.6ti:^ 


c 4+U.V 


115.6 


2.2 


15.6 


4+^2 


llt^ 


0+8 

»_4 


1 


2 


66+°-34 

U.00_o 22 




4.3+1? 


99.8 


8.1 


17.8 


13lf 


54 


5+9 

3-3 




3 


65+° i3 


4+1.6 


q n+0.9 
"^■"-0.8 


98.9 


8.0 


15.2 




5+3 

3-2 


34 




4 


n 09+0-35 
'-'•^^-0.26 


n+1.4 
"J-'-'-O.O 


2 fi+1-2 

'^■"-0.7 


103.8 


9.9 


14.3 


3041 


4+3 


2 O+i "^ 




5 


^•' '-0.31 


9 q+1.9 
^■"-•-0.9 


1 t: + l-7 
^■^-1.5 


87.2 


9.1 


12.3 


qq + 56 
-J '-'-22 


a4 


1 c+2.5 
^■3-0.6 


1 


6 


1 .0+0 AO 
l-^"-0.30 


q 4+1-5 


9 q+0.8 
^•'-'-0.8 


107.4 


11.6 


15.0 


35^?° 


a4 


9 9 + 1.5 
^■^-0.6 


1 



b = 500 km s-^ (fcrkin = 72 keV) 



1 


11+0-03 
^•-'^-'^-0.03 


1 14+0-13 
-1-^^-0.12 


5 0+'^-*^ 

3-"-0.6 


115.6 


1.9 


8.5 


10+^ 


5.6 


2 


29+° °^ 

"•"^^-0.06 


1 97+0.24 
^•^ ' -0.20 


4.0+1? 


99.8 


4.7 


9.7 


23+12 


3.3 


3 


sn+0-06 

"•"^^-0.03 


1 n+0-17 

l-l^ -1-0.15 


2 6+1-1 


98.7 


4.9 


8.3 


27+1° 


3.8 


4 


36+0-08 

"■«J"-0.06 


1 nc+O.lS 
1-U3_Q 


9 9+1.3 
^■^-0.7 


103.7 


5.4 


8.0 


34ti^ 


3.5 


5 


•54+0-10 


95+0-23 
u.yo_Q 20 


1-1-1.5 


86.3 


5.4 


7.7 


35t?3 


3.6 


6 


n 41 +0-08 

^■^-1-0.07 


1-1^-0.15 


1 6+1-° 

l-"-0.3 


108.3 


6.5 


8.5 


374t 


4.0 



-1.2 
+1.1 
-0.9 
+0.8 
-0.7 
+1.0 
-0.8 
+1.6 
-1.1 
+1.1 
-0.8 



-1 q+0.6 

l-'-'-0.3 

u.y_g 2 

0.8 ±0.2 

7+°-2 
^■' -0.1 

0.6+°-2 



0.7: 



0.1 
0.2 
0.1 



2000 km s-^ (fcrkin = 1200 keV) 



1 


11+0-01 

'-'•11-0.03 


0.89t°:°^ 


5-04. 


123.7 


1.9 


7.2 


121^ 


4.7 


2 


26+0-05 

"•^"-0.04 


0.97t°:i3 


4.011:1 


95.0 


4.1 


8.5 


271? 


3.0 


3 


o.26t°:°l 


n 87+0.09 

L'-^'^-o.oo 


2 5+1-° 


101.1 


4.5 


7.1 


304 


3.6 


4 


0.29t°:°l 


0.84t°:°^ 


2 2+1-1 


100.2 


4.9 


7.5 


3511° 


3.5 


5 


o.28t°:^I 


0.80i°:}3 


-1 -1+0.9 
1-1-1.6 


82.6 


4.4 


6.9 


3511? 


3.7 


6 


n q9+0.04 
'-'•''^-0.03 


o.87t°:°« 


1 6+1-1 

l-"-0.8 


110.4 


5.1 


6.9 


37+? 


4.0 



-0.5 
+0.7 
-0.5 
+0.7 
-0.6 
+0.6 
-0.6 
+ 1.0 
-0.9 
+0.6 
-0.6 



0.71°:? 

0.6±0.1 
0.53 ±0.06 

0.5±0.1 
0.55 ±0.06 







Notes. Errors represents 90% confidence limit of statistical errors. See \ 
width of He-like and H-like iron K/3 under the best-fit KABS model, t 
values, t ^-parameters and total column density are calculated based on ( 
are rejected based on overestimation of K/3 equivalent width (denoted as 1) or large A^tot as > 1.6 X 10^'' cm~^ (denoted as 
2). II Plus is defined as blue shift. * Degrees of freedom. 



5 2.3 for systematic errors. * Predicted equivalent 
Ratio of N-pe XXV to -^Fe XXVI from the best-fit 
I equations (1) and (2) derived in §4. § The models 
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Fig. 8. Relative ionization population, /iom of iron (a) and the ratio of He-like to H-like iron (b), shown as a function of the 
ionization parameter ^. Photoionization-recombination oquiUbrium is assumed in an X-ray irradiated matter with a small optical 
depth. The irradiating X-ray spectral shape is assumed to be a DISKBB of fcTin = 1.4 keV (solid line) and fcTjj, = 1.2 keV (dashed 
line) . 



The ratio can be well approximated by a power law, as 
■^Fe XXV /Nfb XXVI ~ 4 X 10^ •<^~^ '^. Similarly, we obtain 
the fraction of H-like iron as /xxvi ~ 1.2 x 10"* • ^^O-^® at 
high values of <^ (roughly ^ > 3 x 10**). Since a value of 
TVtot is obtained as iVtot ■ ^Fc • /xxvi = ^Fo xxvi, where 
j4fc = 3.3 X 10"'"' is the solar iron abundance, it can be es- 
timated as Ntot = -^Fe xxvi/ (^Fo • /xxvi)- By using these 
relations, we derive ^ and iVtot for the range of b given in 
Table 3. The derived values of ^, TVtot and Nxxv/Nxxvi 
are also summarized in the same table. In general, all 
these quantities are anti-correlated with b, and at low ve- 
locity dispersions of 6 = 42 and 100 km s~^, A^tot exceeds 
1.6 X 10^^ cm~^ i.e., is optically thick to electron scat- 
tering. This case is indicated as "2" in the last column 
("Note") of the table. Thus, the case of 6 < 100 km s"^ is 
ruled out by both the K/3 absorption structures and the 
condition of optically thin scattering. We can thus con- 
clude that there must be higher velocities present than 
expected simply from the temperature of the material, 
pointing to the importance of either a structured veloc- 
ity field as expected in an accelerating outflow, and/or 
turbulence. 

4-.2. Time evolution of the absorber 

We repeat the above analysis for all the datasets. Since 
the value of fcTin changes as the source decline, we cal- 
culate XSTAR models Tin spanning the observed range 
of 1.2, 1.3, and 1.4 keV, and then use this to inter- 
polate the ion populations for each observed tempera- 
ture. The dashed lines in figure 8 show the results with 
kTin = 1.2 keV, in comparison with the solid lines which 
assumes kTin = 1.4 keV. We approximated the behavior 
of the ratio of He-like to H-like ion columns and /xxvi by 
following two formulae, 

= 4xlO^-(fcrin/1.4)-5-57.^-l-4 ^(1) 



log(/xxvi) = -2 + 1.69 • exp 



-(log(?/eo))^ 



(2) 



Fo XXv/A^Fo XXVI = 



where /xxvi peaks at log(Co) = 4.13 • (/cTi„/1.4)-0-53. 
Using these formulae, we converted the observed values 
of iVxxv/A^xxvi and Axxvi into A^tot and ^, as given in 
Table 3. For all the data sets, 6 = 42 km s'^ (fcTkin = 
0.5 keV) and b = 100 km are again ruled out because 
A^tot becomes optically thick, as well as generally over- 
predicting the iron K/J absorption (§ 3.3). 

The inferred time histories of the values of ^ and A'tot 
are shown in figure 7 assuming b = 500 km s~^, together 
with upper and lower bounds corresponding to 5 = 200 
and 2000 km s~^. Assuming that the velocity dispersion 
stays constant at b — 500 km then there is a marginal 
decrease in A'tot as the luminosity decreases, from ~ 1 x 
10^'^ cm~^ in the first observation to ~ 7 x 10^^ cm^^ in 
the final dataset, as well as a marginal decrease in ^ from 
- 6 X 10'' to - 4 X 10''. 

4-3. Location of absorbers 

We can roughly estimate n and r of the absorbing 
plasma, by using the definitions of Atot = Jfi^ n • dr and 
^ = L/nr"^ where Rq is the launching radius of the ab- 
sorber. In this section, we use b = 500 km s^'^ to estimate 
^ and A'tot, and choose to replace luminosity with flux, 
so as to make the dependence on the poorly constrained 
source distance explicit. We parameterize this in terms of 
I?io, denoting distance in units of 10 kpc, and note there 
is no uncertainty in ^ parameter associated with the incli- 
nation, since the absorber is along the same line of sight 
as we are. 

The description of A^tot can be simplified as A'tot = n/S.R 
with a mean number density n and a characteristic thick- 
ness of the absorber Ai?. Using ^ = L /nR^ , the average 
distance of the absorber, R, is thus estimated as i? ~ 3 x 
101° • (AR/R) ■ Dio^ and ~ 4 X lO^" • (AR/R) ■ cm, for 
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o N , 
• N, 



T — \ — ' — ' — ' — ' — \ — ' — n ionization state throughout the matter. We thus assume 

^ that n drops as to keep ^ constant at different dis- 

T ; tances. This is also physically realistic for a geometrically 

; thick absorber, since mass conservation requires this de- 

I I * pendency for a constant outflow velocity with constant 
opening angle. In this case, A^tot and ^ are described as 
-^tot = noRa and ^ = L/noRo^, respectively, where uq is 

^ the launching density at Rq. The launching distance is 
thus estimated as Rq ~ (3-4) x 10"'^° • Dio^ cm throughout 

J the observations, with the associated density of no ^ (5- 

1 2) X 10^^ • D^Q cm""^. These values arc in good agree- 

^ ' ' ' ' ^ ' ' ' ' ^ ' ' ' I ' - ment with the previous estimation of the mean distance 

^ j I I ^ J i? as expected as averaging the density profile gives an 

f ; effective width of AR/Rq ~ 0.7. There are no significant 

H — ' — ' — ' — ' — \ — ' — ' — ' — ' — \ — ' — ' — ' — ' — \ — ' — ^ changes over the six observations, but the data are con- 
sistent with a small increase in Rq and/or a small de- 

+ - - + crease in uq as the luminosity decreases. There are some 

\ , , , , \ , , , , \ , , , , \ , uncertainties on these numbers due to the statistical un- 

; certainties in A'tot and ^, but the larger source of uncer- 

~ ^ ^ ~ ~ _ _ - tainty is systematic from the assumed velocity dispersion 

^ ~ ' ' of 6 = 500 km s^^. The lower limit on b gives an upper 

1_ i ^ i ; limit on A^tot and ^, and thus results in smaller launching 

H — ' — ' — ' — ' — \ — ' — ' — ' — ' — \ — ' — ' — ' — ' — \ — ' — ^ radius with larger density. Considering our range in b of 

, " 200-2000 km s~^, the allowed range of i?o is estimated to 

be (0.2-4) xlQi" • Dio"^ cm and (l-5)xl0i° • Dio^ cm for 

* 1 the first and the last observations, respectively. 

190 200 210 220 

Days since 2005 Aug. 13 5. The nature of the wind 



Fig. 7. Time histories of absorption line parameters based 
on the Voigt profile fitting for 6 = 500 km s~^. The top panel 
shows column density of H-like (open circle) and He-like iron 
(filled circle) in unit of 10^* cm~^, the second panel shows 
^Fe xxv/^Fe XXVI i the third panel shows the line blue 
shift, 2, in unit of 10~^. Systematic uncertainty of ~ 20 eV 
for the absolute energy is interpreted as the uncertainty of 
2 as Az ~ 3 X 10~^ at 7 keV, which is also shown in the 
third panel. Based on the ratio, ^-parameter and the total 
hydrogen column A'^tot are calculated and are shown in the 
next two panels in units of 10* and 10^^ cm~^, respectively. 
The best estimation of ^ and A'^tot for b = 200 km s~^ (dashed 
line) and 2000 km s~^ (dash-dot line) are also shown in the 
same panels. 

the first and the last observations, respectively. The den- 
sities are also estimated as 5 x 10"'^^ • (AR/R)^^ ■ Dio~^ 
and ~ 2 X lO^^ . (^^Ji/Jiy^ . d^-'^ atoms cm . Assuming 
Ai?/i?~ 1 by following the previous research (e.g., Kotani 
et al. 2000; Ueda et al. 2001; 2004), the values of R and n 
are estimated to be on the order of i? ~ 10^" • Diq^ cm and 
fi~ 10^^-Dio~^ atoms cm~^. They depend on the plasma 
thickness AR/R, and the longer distance and higher den- 
sity are realized on the the condition of smaller AR/R. 
However, since this structure is continuous, it seems un- 
likely to have too small thickness (e.g., Ai?/i?< 0.1). Thus 
the range of the distance can be securely estimated as 
igs-io . This value is well consistent with those 

estimated from the other binaries which showed absorp- 
tion line structures. 

Because of lack of any evidence for low ionization mate- 
rial, it is reasonable and convenient to assume a constant 



The estimated launch radii are much smaller than the 
distance traveled by the absorbing material over a sin- 
gle 20 ks observation from the measured blue shift of 
^ 1000 km s~^, 10^^ • DIq cm. Thus the persistence of 
the wind over the two-month monitoring period, where 
the source was always in the high/soft state, shows that 
the material is not associated with a one-off event such 
as a state transition. Neither can it be connected to the 
jet since the radio emission is strongly suppressed in the 
high/soft state (Fender et al. 1999). Instead it must be 
produced by the disk in a continuous or quasi-continuous 
manner. 

Here we use the estimated physical parameters of the 
absorber to constrain the driving mechanism for this wind. 
One of the candidate mechanisms is radiation pressure on 
electrons, which becomes dynamically important as L ap- 
proaches iEdd, reducing the effective gravity by a factor 

(1 — L/LEdd)- This can be made much more efficient if 
the cross-section for interaction between the matter and 
radiation is enhanced by line opacity. There are multi- 
ple line transitions in the UV region of the spectrum, so 
UV emitting disks can drive a powerful wind at luminosi- 
ties far below Eddington. Such line driven disk winds are 
seen in CV's (Pereyra et al. 2000) and are probably also 
responsible for the broad absorption line (BAL) outflows 
seen in AGN (Proga. Stone, Kallman 2000). However, 
the disk temperature for black hole binaries is in the X- 
ray regime, so line driving is probably unimportant here 
(Proga & Kallman 2002). We do indeed see lines, but 
the momentum absorbed in these iron transitions is very 
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small, completely insufficient to drive the wind. 

Another type of outflow from a disk is a thermally 
driven wind (Begclman et al. 1983). Here again the cen- 
tral illumination is important, but the process is less di- 
rect. The illumination heats the upper layers of the disk 
to a temperature of order the Compton temperature, Tc- 
This will expand due to the pressure gradient, and at large 
enough radii the thermal energy driving the expansion is 
larger than the binding energy, leading to a wind from 
the outer disk. Simple estimates of the launching radius 
of this wind give R = 10" • (M/Mq) ■ {Tq/IO^ K)-i - lO^^ 
cm (Begclman ct al. 1983). This is a much larger distance 
than inferred for our wind. However, this is an overes- 
timate, and a more careful analysis shows that thermal 
winds are launched at a radius a factor 5-10 smaller than 
this (Begelman et al. 1983; Woods et al. 1996) . Secondly, 
the combination of thermal driving and radiation pressure 
from moderately sub-Eddington illumination can be much 
more effective at producing a wind. The hydrodynamic 
calculation of Proga & Kallman (2002) shows that ther- 
mal driving with L/L^dd ~ 0.6 gives a powerful wind from 
well within the expected thermal driving radius. Their 
simulation is not exactly matched to the situation here, 
as they model a neutron star rather than a black hole. 
Nonetheless, their wind is launched from a distance of 
around 1000i?s, has a column between 10^2-23 qyr~'^ and 
ionization of ^~ 10^ for inclinations larger than 70°. This 
is very close to our inferred parameters for the wind seen 
here in 4U 1630 — 472. The simulation even has an outflow 
velocity of ~ 1000 km as currently indicated by the 
data. 

The last type of outflow is a magnetically driven wind. 
These arc much harder to quantitatively study as the 
magnetic field configuration is not known, yet they are 
almost certainly present at some level as the underly- 
ing angular momentum transport is known to be due to 
magnetic fields (see e.g. Balbus & Hawley 2002). Winds 
(and jets) are clearly present in magnetohydrodynamical 
(MHD) simulations which include these magnetic stresses 
self consistently. These generically show that the mass 
loss is stochastic, with large fluctuations both spatially 
and temporally, but that the time averaged properties are 
well defined, so this magnetic wind is quasi-continuous 
(e.g., Hawley & Krolik 2001; Machida, Nakamura, & 
Matsumoto 2004). These calculations are still in their 
infancy, especially for describing the properties of a ge- 
ometrically thin disk as appropriate here, so cannot yet 
be used for quantatiative predictions of the wind in the 
high/soft state. Instead, an approximation to the proper- 
ties of the self-consistent magnetic wind from an accretion 
disk can be made by imposing an external field geometry 
(Proga 2000; 2003). The mass loss rates depend on this 
field configuration, but in general these allow steady, pow- 
erful winds to be launched from any radius. This lack of 
diagnostic power means that they can be reliably identi- 
fied only when all other potential mechanisms for a wind 
are ruled out. 

While the physical parameters required for the absorb- 
ing material seen in 4U 1630 — 472 are somewhat more 



extreme than produced by the simple analytic estimates 
of a thermally driven wind, they seem to be potentially 
within the range of composite thermal/radiation pres- 
sure driven winds if L/L-^dd > 0.1. Unfortunately there 
are no firm estimates for the distance and black hole 
mass for 4U 1630 — 472, and thus the Eddington ratio 
(L/LEdd ^ 0.2-0.3 • a) is hard to constrain. Neither are 
there any directly applicable simulations of these winds, 
and how their properties might change with L/L^m de- 
creasing by a factor 1.7. Hence there is no unambigu- 
ous requirement for magnetic driving from these data. 
Nonetheless, this wind is clearly very similar in properties 
(though with a slightly higher ionization state) to that 
that seen in a similar high/soft state from another black 
hole, GRO J1655-40, where Miller et al. (2006) inferred 
a very small launching radius from the Chandra grating 
data requiring a magnetic wind. Given the continuity of 
properties seen between the wind in 4U 1630 — 472 and 
GRO J1655 — 40 then it seems likely that if the wind in 
GRO J1655 — 40 is magnetically driven, then this process 
is important also in the 4U 1630 — 472 absorber examined 
here. 
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